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1 Introduction and summary 



One of the hallmark features of open bosonic string field theory is the existence of a 
tachyon vacuum state, around which there are no perturbative excitations. The perturba- 
tive vacuum may describe any consistent D-brane configuration, depending on the choice 
of boundary conformal field theory (BCFT). As explained by Sen [2], open bosonic string 
field theory (OSFT) built upon arbitrary BCFT possesses always another vacuum which 
corresponds to a state with no D-branes and hence no open string dynamics. Many other 
solutions have been constructed since then, either numerically or analytically. Some of 
these solutions describe lower dimensional D-brane configurations. From the viewpoint of 
the theory expanded around the tachyon vacuum, these solutions correspond to various 
D-branes popping out of the vacuum, but up to date all solutions found have less energy 
than the original D-brane. On one hand it might not appear so surprising, as one expects 
the tachyon condensation to drive the D-brane system to a state with smaller energy. 
On the other hand, since the final no-D-brane vacuum state is believed to be unique, 
and string field theory formulated around this vacuum does have solutions with positive 
energy, it is clear that the apparent impossibility to go higher in the energy is akin to the 
insufficiency of a particular coordinate system to describe the whole geometry. 

The question is thus, how big is the space of string fields formulated around a given 
reference BCFT and whether solutions with positive energy with respect to the pertur- 
bative vacuum do exist. This issue has been partially numerically studied in the past by 
Taylor, Ellwood and also by the second author, using level truncation method, but no 
conclusive evidence for the existence of such solutions was found. An analytic solution 
was proposed by Ellwood and the second author, but attempts to compute its energy in 
level truncation yielded a number that was off by roughly a factor of minus twelve. 

In this paper we are going to study a class of universal solutions of the OSFT equations 
of motion in the form 

KB 

m = Fc-^-^^cF, (1.1) 

where F = F{K), and K, B, c are, by now, well known string fields. These solutions are 
all universal, in the sense that their form does not depend on the detail of the BCFT. 
They do exist on any D-brane configuration. Analytic solutions of this form have been 
studied in [31 lU El El [7] but always with additional assumptions on the function F which 
more or less guaranteed that the solution was tachyon vacuum. Our proposal is to adopt 
the least possible assumptions on F and compute the energj@ and the Ellwood's gauge 

^The energy for a class of tachyon vacuum solutions of this form with F{0) — 1 and F'{0) > was 
first computed by Erler in [7]. His computations, following quite closely ^ and being perturbative in F, 
are unfortunately not general enough for present purposes. 
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invariant to deduce the physical meaning of such a solution. 

As this paper is rather technical, let us summarize our main assumptions and results. 
It turns that the appropriate conditions can be more conveniently stated in terms of a 
complex function 



The most elementary, but important condition is that G and 1/G have inverse Laplace 
transform which has non-vanishing support only on the positive real axis. This means 
that we can introduce the function G{K) of a string field K 



as a superposition of non-negative width wedge states e~ and hence the star product 
of such states can in principle be well definedjj This condition can be restated as the re- 
quirement that both G and 1/G are holomorphic and that their absolute value is bounded 
by a polynomial in the half planes Kez > e for any e > 0. As we shall see this condition 
is not very strong and will allow occasionally for some divergences, and in particular it 
would not guarantee that the energy computed from the action and from the Ellwood 
invariant would agree. 

Stronger conditions which simplify some computations and give unique answer for the 
energy can be summarized as follows: 

i) G and 1/G are holomorphic in Rez > except at 2; = 0. 

ii) G is meromorphic at z = 0. 

iii) G is holomorphic at z = 00 and has a limit G(oo) = 1. 

With the help of these assumptions we can evaluate easily the action and we find for the 



where the closed contour C encircles all singularities and branch cuts in the Re^; < 
half plane. It does not wind around the origin however. Had some of our assumptions 
been violated, e.g. the integrand had branch cuts that extend to the boundary of the 
half-plane, the energy could still be computed, but appropriate limits must be taken (if 
they exist) and additional terms might appear, as we shall discuss in Sec. [31 

**More detailed discussion can be found in [8]. 






(1.3) 



energy E = 1 ( ^, Q^^ ) of the solution (dUD 




(1.4) 
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In our second, independent computation we evaluate the energy using the EUwood 
invariant 

£ = (/kavWI*) = -^i.m^f|. (1.5) 

where V^r(^) is the properly normalized zero momentum graviton vertex operator inserted 
at the string midpoint i and |/) is the identity string field. How can the two expressions 
(11. 4p and (11. 5p be compatible? Our stronger condition ii) guarantees that the function 
G' /G is meromorphic in z = and hence (II. 5p can be written as a tiny contour integral 
around the origin. At the same time G' / G is holomorphic at infinity and this allows us to 
deform one contour into the other to prove that both expressions are actually the same! 
One notable example of a family of functions obeying the stronger conditions is 



Gn{z)={j^] , (1.6) 



for which the energy computed either way is 



n 



B = -^. (1.7) 

For n = 1 this solution represents the tachyon vacuum solution [9], while n = gives 
the perturbative vacuum \1/ = 0. Negative values of n describe states with energy higher 
than the perturbative vacuum. In fact, we conjecture that they describe configurations 
of multiple D-branes. Positive values of n would describe "ghost" branes, objects with 
negative tension. Such objects are not expected to arise in bosonic string and we show 
that they are indeed divergent in Fock state expansion. 

This paper is organized as follows: In section |2] we develop some tools that enable 
us to compute the energy of our solutions quite efficiently, and we do it in a number of 
ways. In the following section we discuss possible anomalies when some of our stronger 
conditions are not met. Section H] is devoted to the study of the solutions from level 
expansion perspective. We will demonstrate that our solutions are well defined in level 
expansion, nevertheless, we will present arguments for the necessity of adding a so called 
non-vanishing phantom term. We then test numerically one specific proposal, suggested 
by several analytic computations, but we do not reach definite conclusions. Section |5] is 
devoted to the discussions of various issues touched upon in the main body of the paper. 

Note added: This is an extended and more detailed version of our conference pro- 
ceedings report [lO] . As the current paper was nearing the completion, two very interesting 
papers appeared: [HI [12] . The first one computes the boundary state for the solutions of 
the type we study, while the second one discusses related results from a broader geometric 
perspective. Both papers report similar difficulties which we believe are due to not yet 
fully understood phantom terms. 
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2 Computations of the energy 



2.1 Useful correlators 

To compute the energy using either kinetic or cubic term of the action it is necessary to 
evaluate number of correlators of the form 

{F,,F2,F,,F,) = {F,{K)cF2{K)cF,{K)cF,{K)cB) , (2.1) 

where Fi{K) are ghost number zero string fields, given by arbitrary functions of K. In this 
paper we restrict our attention to the so called geometric string fields in the terminology 
of [8]. This condition means that Fi can be written as formal Laplace transforms 

F,{K)= / /.(«)e-"^, (2.2) 



of arbitrary distributions fi with support in [0,oo). One can impose various additional 
conditions on the distributions as discussed in |8], but in this paper we will formulate 
them in terms of the properties of the functions F. 

The advantage of such a restriction is twofold: on one hand it gives us a nice geometric 
picture of the string field Fi[K) being superpositions of wedge states e"""^ (note that a has 
to be non-negative), and on the other hand it offers means of computing the correlators 
through the formula [Sj HJ [7] 

(e-"i^ce""2^ce-"3-^ce""4-^cS) = (c(ai)c(ai + a2)c(ai + 02 + a3)c(ai + 02 + as + "4)5)^^ 



4v]-3 



. 27ra2 / X . 27r(a2 + as) 

04 sm (as + 04) sm 



s s 



2Tra4 271(0:3 + 04) 

+a2 sm (^2 + as) sm 

s s 

. 27r(a2 + ag + 04) . , , . . 2Tia^ 

+03 sm h (02 + 03 + 0:4) sm 

s s 

which expresses the left hand side as a correlator of four boundary c-ghost insertions 
and one line-integral 6-ghost insertion on a unit disk presented as a semi-infinite cylinder 
Cs of width s, with the midpoint mapped to infinity. To shorten the notation, we have 
introduced s = Ylt=i ^i- explicit formula for the correlator (12. ip can be obtained from 
the quadruple integral 

/•oo roo /"OO roo 

/ / / / daida2rfa3rfa4/i(ai)/2(a2)/s(as)/4(a4) (e-"^-^ce-"2^ce-"^^ce-"*^c5) 
Jo Jo Jo Jo 

(2.3) 
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by a simple trick. Let us insert into the integral an identity in the form 



/ (is 5 s — > ai ] 



,+ioo ^. 



J —too 



ds / e'"e-"^-i"% (2.4) 

2Tci 



which allows us to treat s as independent of the other integration variables at. The 
second equality is just the ordinary Fourier representation of the delta function with the 
i absorbed in the integration variable, so the contour runs along the imaginary axis. 
The integrals over can be easily performed and reexpressed in terms of the original 
functions Fi{z) 



F,A{F2F;)Fi - F.F^FsAF, + F,F^A{F,F,) + F.F^AiF^F^) 



-F,A{F2F^F,) - Fi{F2AFsF^y 



(2.5) 



where for convenience we have omitted common argument z and also introduced an op- 
erator As defined as 

(A.F)W = F(.-^)-f(.+ ^), (2.6) 

When no confusion arises we omit the subscript s. Let us now list some of the properties 
of this quadrilinear correlator: 

(Fi,l,F3,F4) = 

(Fi,F2,l,F4) = 

(Fi,F2,F3,l) = 

{F,,K,K,F,) = 

{FuF2,K,K) = 

{K,F2,K,F,) = (2.7) 

The first three equations are true because = and the second two express cKcKc = 0. 
All of the five equations are obeyed by the integrand of (12. 5p itselfjf] The sixth equation 
can be written as 

{QBiBc)F2QBic)F,) = (2.8) 



^The fact that the integrand itself manifestly obeys these five conditions is the reason why we prefer 
this correlator to the one introduced in 1101. 
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and hence should be true by the basic axiom (Qb(...)) = 0. The square-bracket part of 
the integrand in fl2.5p for the choice of Fi = F-^ = K is equal to 



z{AF2Fi+F2AFi-A{F2Fi))-z^AF2F!,+zA{zF2)F!,-z^F^AFi+zF^A{zFi)-z{F2AzFiy. 

(2.9) 

The discrete derivative obeys a sort of deformed Leibniz rule 



s J \ s J \ sy \ s J 
= A2s{fo.,g). (2.10) 

With the help of this identity f l2.9p can be rewritten as 

A2s{F2 o, zF^) - A2s{z o, F2F4) - A2.(F2 o, z^f:^) + A2s{zF2 o, zFi) + 
+ -^{zd, - sds)i-sF2AlF,) . (2.11) 

Now let us consider the ^-integral in f l2.5p . If F2 ■ F4 is at most 0{z^) at infinity, then 
the integrand is no worse than 0{z~^) and so we can turn the 2;-integral into the integral 
over a closed contour Cg by adding a noncontributing arch at infinity in the left half plane 
Re 2; < 0. The subscript s reminds us, that when we try to shrink the closed contour to 
a finite one, the minimal dimensions of such a contour will typically depend on s. The 
discrete derivative A2S vanishes under the integral sign 

/ -^e-A2.(/io,/2) = 0, (2.12) 

for any pair of possibly s-dependent functions fi^2{z), provided that the closed contour 
Cg is sufficiently large so that it encircles all the singularities of fi{z) and fi{z ± 2iTi/s) 
in the left half plane Re ^ < 0. We thus arrive at 

{K,F2,K,F,) = ds £ — e- — (^9, - sdg){-sF2AlF,) . (2.13) 

The operator {zdz — sdg) commutes with the exponential e'^^ and we find that the whole 
expression reduces to possible surface terms 

{K, F2, K, F4) = f lim - lim) / dz e''-^s^F2Al,F4 . (2.14) 

Vs-s>oo s-s>0/ J^^ OTT'^l 

The surface term at s = vanishes if both F2 and F4 are at most 0{z) at infinity and the 
one at s = 00 vanishes if F2d'^F^ does not have poles on the imaginary axis. 
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2.2 Kinetic term 



Let us now describe our computation of the energy using the kinetic term. By a simple 
manipulation we are led to consider the sum of the following four correlators 



vl>,Qv]>) 



K K 



(2.15) 



^,{1-G),K,k\ + Ik,{1-G),K,^ ). 



The third term vanishes by (12. 7p . The fourth term can be omitted as well using the 
stronger assumptions on G, but for the sake of generality we shall keep it. By applying 
the s-z trick (12. 4p . the kinetic term can be expressed as 







ds i - — - 



dz 



27ri Svr^i 



Aiz^G'^ 



2z'G'A 



G \G 
which can be further simplified using the Leibniz rule (I2.10p 

dz e' 



(2.16) 



ds 



2-iTi Sn^i 



2A'Kisz^% - 3(za, - sds) (s^z^^^ 
G \ G 



(2.17) 



+ 2s^A2s ( z o 



z'G 
~G 



s'^Ao.s ( zG o 



z'G 
~G^ 



If G is meromorphic at z = cxo, that is G can be written in the form G = f Xl^o '^n^ "' 
then the square bracket part of the integrand decays as 0(1/^^) or faster. With this 
condition one can now make the integral along the imaginary axis into a sufficiently large 
closed contour Cg by adding a non-contributing arch at infinity in the left half plane 
Rez < 0, 



1 fd. / ^ 

Jo Jcs 27^2 



sz^^-^Xzd.-sd.) (s-z^ 
G 8m \ G 



(2.18) 



where we have dropped the second line in (I2.17P by using (I2.12p . 

The integration over the remaining two terms can be performed separately. In the 
first term the Gs contour can be chosen to be s independent, as the singularities do not 
depend on s. The double integral then makes perfect sense, so the s integration can be 
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performed firsll§ and we find 

dz G'(z) 



7r2 27ri G{z) ' 



(2.19) 



We will see that this value is consistent with the one coming from the Ellwood's gauge 
invariant observable if G{z) satisfies the conditions i)-iii). 

Therefore the second term of fl2.18p . if nonzero, is an anomaly. In this term, we can 
move the zdz — sdg operator outside the exponential factor. Integrating by parts, only a 
single surface term can contribute 

Now, since the contour has been chosen generously large, its infinitesimal variation does 
not change the integral and hence one can move the derivative outside the integral, giving 
rise to the two surface terms 

3 / . \ I dz izs^ ^izGMz) 

— (lim-lim)i e'^ , ' . 2.21 

These surface terms could give anomalous contributions so that the value of the kinetic 
term would be inconsistent with the gauge invariant observable. Under certain conditions 
and/or a prescription given in Sec. El these contributions vanish and we find 



, ^ ^ ^ , 3 f dz G'iz) 



which establishes the result (II. 4p . Note that the integral is the well known topological 
invariant on the space of meromorphic functions, which counts with multiplicity the num- 
ber of zeros minus the number of poles inside the contour G. In our case, however, we 
allow the function to have arbitrary, even essential singularities inside the contour, but 
nothing outside except for a possible pole or zero at the origin. 

One of the crucial assumptions is that both G{z) and z/G{z) are Laplace transforms 
(corresponding to geometric string fields), and hence G does not have any poles or zeros 
in the half-plane Rez > 0. The little bit stronger conditions i)-iii) allow us to shrink the 
G contour around infinity, picking up only a possible contribution from the origin. 

It dzG'(z) 



^ If we did not encircle the singularity oi G' / G{z) at zero from the left, the double integral would have 
depended on the integration order. 
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where Co is now a contour around the origin. There is no contribution from infinity if it 
is a regular point. For G{z) ~ near the origin we get 



n 



(2.24) 



E 



Thus the integer n + 1 can be interpreted as the total number of D-branes on top of each 
other, including the original D-brane. 

At first sight there seem to be no reason why poles of G at the origin would be allowed 
while zeros not. We will see in Sec. 14. II that for n < —4 the solutions are singular in level 
expansion. While for the solutions with n > 3 some of the coefficients are also singular, 
we believe this is a milder singularity which will be canceled by a phantom term discussed 
in Sec. 14.21 The solutions with negative value of n have a tension which is too negative 
and should correspond to the unphysical states with negative number of D-branes. We 
will call such states as the ghost branes, although this term has been already used in the 
literature in a rather distinct context. 

2.3 Energy from the cubic term 

The energy of a string field theory solution can also be computed using the cubic term 



Often, this is a useful check, since the solution might not obey the equation of motion 
automatically when contracted with itself. That one gets the correct answer for the 
tachyon vacuum in the Bq gauge was verified successfully in [U |5] . 

Using the explicit form of the solution (11. ip and reducing the number of B insertions 
by (anti) commutation we arrive at 



Now using mere linearity of the four bracket and the identity KG/G = K, which is 
trivially suggested by adopting the s-z trick (i.e. it holds at the level of the integrand 
(12. 5p ). we get exactly the same term as from the kinetic term (12.150 . Hence there is no 
new computation to be done. 



E = — (^,^*^). 



(2.25) 



(^,^*^) = / 1 




(l-G)^ 
(2.26) 
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2.4 Ellwood's gauge invariant observable 

We shall evaluate the energy through the gauge invariant observablei 



{I\V{m), (2.27) 

discovered by Hashimoto and Itzhaki [H] and independently by Gaiotto, Rastelli, Sen 
and Zwiebach [TBj. They depend on an on-shell closed string vertex operator V = ccV"^ 
inserted at the string midpoint i in the upper-half-plane coordinates. The state (J| is 
the identity string field and it also plays the role of the Witten's integration. EUwood 
conjectured [13] that the gauge invariant observables compute the difference of the one- 
point functions of the closed string on the unit disk between the trivial vacuum and the 
one described by the solution \Ef: 

(/|V(z)|^) = Af^V"") - Af'^V^) . (2.28) 

Here ^^^'^'^(l^™) is the one-point function of the matter part of the closed string vertex 
operator on the disk in the vacuum \E'. For the tachyon condensation solution, Ellwood 
confirmed that the one-point function is zero. For the A^-brane solution the one-point 
function is expected to be N times the one in the trivial vacuum: ^^''^'^(V^™) = NAq^^^^V^). 

Now let us compute the gauge invariant observable for the solution (11. ip . Since K 
commutes with the operator inserted at the midpoint, the invariant equals 

(vFcBFcF^ = (VcBFcF'^^ , (2.29) 

where F = K/{1 — F"^) and V = V(i)|/). When we express and F as the Laplace 
transform and denote the inverse Laplace transforms as / and / respectively, it becomes 

roc /»oo 

da d(3f{a)f{(3){VcBe~''^ce-^^). (2.30) 
Jo Jo 

A very similar correlator has been computed by Ellwood [13], so with the help of a simple 
reparametrization we find 

^OO poo _ r)- _ 

/ da dPf{a)f{P)-P{VnvT= ^^F{z)d^F'iw)Af^\vn 

Jo Jo TT z,w^O 

= -(piS^f|)-4»"(''"'' P.31) 

''We can evaluate the boundary state following [16], but for simplicity we will evaluate the gauge 
invariant observable only. 
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where ( ■ )uhp'''^ matter correlator on the upper half plane (UHP) and G = l — F"^. In 

?2l 



general d^F (z) or F{z) may diverge in the z — )■ limit, but their product is much better 
behaved. There is a natural regularization which makes this manifest. Imagine replacing 
string field K with K + e in the solution (11. ip . where e is a small but positive number 
(multiplied by the identity string field). The only change to (I2.30p is the appearance of 
an extra factor of e"^""*"^)^ under the integral sign. After evaluating the correlator we get 
F{e)dF'^{e) and in the limit we justify the final expression in (I2.3ip . 
Substituting the result (I2.3ip into the Ellwood's relation we get 

/ C'(z)\ 



At'-iVn =[l- li^o^^ J X'^'iVn ■ (2.32) 

Now if G{z) behaves as around the origin, 

Af^V"^) = {n + l)Af'^{V'^). (2.33) 

Thus the solution \E' can be expected to describe the (n + l)-brane. In particular, one can 
compute the energy from the EUwood invariant by using appropriately normalized zero 
momentum graviton vertex operator 



27r2 7To G{z) 



E = {I\ccVgr{i)\'^) = -^limzy:^. (2.34) 



For a function G meromorphic at the origin, or more generally one for which G'/G is 
meromorphic, the limit can be replaced by a contour integral 

(2.35) 

where Co is small contour encircling the origin. As we have explained in the introduc- 
tion, this result is consistent with our energy computation from the kinetic term, if the 
conditions i)-iii) hold. 

Let us now show how the result (I2.32p can be obtained directly in the contour form 
using the s-z trick. This exercise will hopefully help the reader to acquire more familiarity 
with the trick itself. By the rules of the trick, the left hand side of (I2.3ip can be written 
as 

Suppose that G = z^g{l/z) with g holomorphic at the origin. Since zG'/G behaves as 
0{z^) around infinity (for r 7^ 0), we cannot simply close the line 2;-integral into a closed 
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contour, but we have to add a delta function 6{s). In total we have 

^d-k(ym) ]^ f^2m G Jc2mG Jc^ 2m G Jco^mG' 

(2.37) 

where Cod is the circle with sufficiently large radius encircling all the poles of G'/G. In 
the last equation, we have assumed the condition i). 

To close the discussion let us see what happens when we combine the s-z trick with 
the K + e regularization. The normalized EUwood invariant fl2.36p becomes 

Assuming for simplicity our condition iii), i.e. G{oo) = 1, the integrand decays rapidly 
in the Re 2; < half-plane and therefore the contour can be closed. The pole (or zero) of 
G at the origin shifted to —e does not contribute to the integral, because of the factor 
z + 6. The contour can be moved left, so that the point —e is just to the right of the 
contour. The resulting double integral is now absolutely convergent, so the s- integral can 
be performed first. We find 

^disk(^m) e^oJc2m z G{z + e) e^o 27fi z - e G{z) Jc^mG ^ ' 

In the second equality we have shifted the variable z to z — e and correspondingly the 
contour by e to the right. In the last equality we took the limit, noting that neither 2; = 
nor z = e contribute, as they lie outside the contour. 

Notice that the invariant comes entirely from the e term in the factor [z + e] in fl2.38p . 
Had we considered only the first z terni§, we would not have been allowed to move the 
contour past — e, and the total contribution would have been zero. This is quite a general 
feature of EUwood invariants that they receive contribution only from a naively vanishing 
piece, a so called phantom term. In Sec. 14.21 we will propose another type of a non- 
vanishing phantom, associated to the residue at 2; = in the s-z trick computation. This 
term does not contribute to the EUwood invariant at finite e. But it does when we set 
£ = strictly at the beginning, in the sense that the contribution from the residue at the 
origin must be subtracted (or not counted) in the solution. 



^This would be actually equivalent to computing the EUwood invariant for a pure gauge solution with 

G{K + e). 
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3 Possible Anomalies 



In the formalism of string field theory one occasionally encounters expressions which are 
either divergent or anomalous and must be treated with due care. One of the least 
understood type of divergences is associated to the so called sliver divergence. The string 
fields e~"'^ describe wedge states of width a, but these do not vanish in the large a limit, 
contrary to what one would naively expect from the large a behavior of e~"^ if K was 
thought to be positive definite. Such anomalies has been recently studied by Erler and 
Maccaferri |T7], following an inspiring work on tachyon lumps by Bonora et al. [18j. In 
our computation of the energy or the Ellwood's invariants, we have repeatedly used the 
Laplace representation of various functions of K, without really worrying whether they 
exist as string fields. This has allowed us to make fast computational progress but leaves 
behind many questions. In our computation of the energy using the kinetic term or the 
Ellwood's invariant we have found two potentially conflicting results. We shall study the 
surface terms fl2.2ip and look under what conditions they vanish. In the last subsection 
we will examine possible anomaly in the equations of motion. 

3.1 Anomaly at s = oo 

Let us look at the s = oo surface term f l2.2ip in the kinetic term f l2.15p 



If the singularities or zeros of G were all at Re 2; < 0, then the contour Cs could be 
deformed, moved away to the left of the imaginary axis and the suppression factor e'*^ 
would guarantee that the s — )• 00 limit is vanishing. For fractional brane^, however, there 
is a branch cut in G going all the way to zero. One can consider for example 



which can be taken to have a branch cut between —1 and for non- integer r. In such a 
case one can suspect that the s = 00 surface term is non- vanishing, which indeed happens, 
as we will now show. Let us write 



assuming that g{z) is holomorphic and nonvanishing at zero. The contour integral in (13. ip 
is well defined for —2 < r < 2. For non- integer r outside this range, the contour integral 

^These singular solutions should not be confused with fractional D-brancs on orbifolds. 




(3.1) 




(3.2) 




(3.3) 
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is still finite, with the prescription of encircling the branch point at z = from the right. 
After a change of variables z ^ z/s the integrand possesses well defined holomorphic 
limit for large s. Closing the integration contour and further changing the integration 
variable to w = we find the anomaly 

^=oo = / T^^^^^—A ((1 - 27rz^)i-^ - (1 + 27iiwf-^) 

= -^(r=^-r)[iFi(2 + r,4,27r2) + iFi(2 + r,4,-27rz)] (3.4) 

This anomaly vanishes only at r = —1,0,1 and at a sequence of non- integer values 
r = ±3.63948, ±8.12955, ±14.2036, .... At these latter values the surface term vanishes, 
but generically we expect the solutions to be still singular. But what about integer values 
|r| > 2? If the anomaly was unavoidable, we would have to conclude that the multiple 
brane solutions (beyond the double-brane solutions) do not exist. For integer values of 
r, one has the option to take the z contour to bypass the singularity at the origin from 
the left0 We will see that the same prescription is needed also to avoid anomaly in 
the equations of motion. This prescription will be interpreted in Sec. 14.21 as an extra 
phantom term visible in level expansion. For non-integer brane with r > 2 this option is 
not possible since one cannot bypass non-integrable singularity on the side of a branch 
cut. Therefore we expect that the equation of motion should not hold if r is fractional 
and so the fractional brane solutions do not exist. 

If G had additional zeros or poles on the imaginary axis, except the origin, and if 
we encircled them from the right, these would give make the surface term (13. ip behave 
as C(s^), i.e. quadratically diwrgentH Therefore, one would conclude that the contour 
Cs should encircle them from the left. However, if we do so, the value of the kinetic 
term becomes incompatible with the EUwood invariants since the zeros and poles on the 
imaginary axis produce anomaly if we deform the contour C into Cq. 

3.2 Anomaly at s = 

Identical computation can be done to study under what conditions does the s = surface 
term contribute. Assuming that the behavior of G around 2; = 00 is of the form 

G{z) = ^rgiljz), (3.5) 

"'^''Note that the integrand of the regular term is regular at zero, so the problem has not arisen in our 
previous discussion. 

^^These divergences cancel in the case of the tachyon vacuum described by the function G(z) = 1 — e^^ 
which has an infinite number of zeros on the imaginary axis. 
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with g holomorphic around zero, one finds exactly the same anomaly (13 ■4p up to a minus 
sign due to the overall sign in front of the s = contribution. 

^=0 = ^(^' - [i^i(2 + r, 4, 2m) + iFi(2 + r, 4, -2m)] . (3.6) 

The anomalies ^4^=0 and As=^ may cancel if the power-law behavior at infinity and at 
zero is the same. To get rid of the anomalies in the equations of motion (discussed in the 
next section), one has to pick a contour bypassing zero on the left and then A<j=oo "will 
automatically vanish. The anomaly As=q will vanish independently due to the condition 
iii), i.e. lim^^oo ^(2;) = 1, which we impose anyway. 

3.3 Anomalies in the equations of motion 

Let us take the simplest candidate multiple brane solutions and regularize it by replacing 
K with K + e 

'^ = ^ G^^K + e) Bc{l-G{K + e)). (3.7) 

This regularized solution can be viewed as a sum of a pure gauge string field with G{K+e), 
plus an additional piece proportional to e. In Sec. 12.41 we have seen that this naively 
vanishing term provides the full contribution to the EUwood invariant. This term is 
reminiscent of the ip^ piece of [3] which have been dubbed a phantom term, since it was 
vanishing in level expansion, but nevertheless it contributed non-trivially to the energy. 
Here, the term linear in e is equally important, and vanishes in level expansion as well. We 
shall call it the e-phantom, to distinguish it from the finite phantom studied in Sec. 14.21 
Thanks to the e-phantom, the equation of motion now reveals a possible anomaly 

As = + ^ * ^ 

and the important question is whether it vanishes as e is sent to zero or not. We shall 
explicitly check its first coefficient in the level expansion. 

To obtain the coefficient tpip2F3 in front of CiCo|0) in a general expression F1CF2CF3 
we need first the elementary formula 

e-'^^ce-^^ce-^^ = ?7i+,_,^+,c g(-a + /3 + 7)) £ (^(-a - /3 + 7)) |0) 

= W7C1C0IO) H , (3.9) 



where 



s + iy f2Y . frcl + 2a\ . /vr 2/3 \ . /vr 1 + 27 
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With the help of the standard s-z trick we readily find 



tFlF2F3 



0) ^r*/_zi^^"(^) [^.(-^)--^.(^+^)-] 

X [F^iz - w) - F^iz + u)] [Fs{z - u)e^ - F^{z + uj)e-^ , (3.11) 



where 



u = ^. (3.12) 

s + l 

Setting Fi = 1, F2 = ^^^^^.^ and F3 = 1 — (^(z + £:), where G is holomorphic in Re z < 
with no zeros or singularities on the imaginary axis except at the origin, we see that the 
2;-integral is completely regularized. It receives contributions from z = —e ± and from 
the points in the interior, i.e. z = — e ± where z^ are zeros or poles of G for which 
necessarily Re z^, < 0. The s integral is then convergent because of the suppression factor 
e*^. In the limit £ — ?■ there might however arise terms 1/e if the non-exponential part 
of the integrand does not vanish for large s. 

To compute such a possible anomaly in the equation of motion we set once again 

G{z) = \g{z), (3.13) 
z^ 

assuming that g{z) is holomorphic and nonvanishing at zero. The integrand of fl3.1ip has 
for r > a pole of order r so to compute the residues one has to perform r — 1 derivatives. 
The leading term in s is obtained when these act on e^^ (2; + e ± ■ After performing 
the residue integral the leading behavior for large s is e~^^ xO{l) giving upon s-integration 
a factor oi 1/e which cancels an overall factor of e. The final answer in e — )■ limit is a 
finite nonzero anomaly 

A=^ (4-i(27r0 + 4-i(-27rz)) ciCo|0) + ■ ■ ■ , (3.14) 

where lI^\z) are generalized Laguerre polynomials, valid for r > —2. Alternative expres- 
sion which makes sense also for non-integer r is 

A = ^r(r + 1) [iFi(2 + r, 3, 27n) + iFi(2 + r, 3, -2m)] CiCo|0) + ■ ■ ■ . (3.15) 
o 

This manifestly vanishes for r = —1,0. For the double brane the coefficient is equal to 
7r/2 for all choices of 5'(2;)0 Since this anomaly is independent of g there is no obvious 
way to cancel it by a clever choice of g. The simplest possibility is to avoid the anomaly 



^^For the 'ghost branes' with r < — 2 there wiU be an additional g-dependent contribution to which 
the first term in = 1 — G contributes. 
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by taking the contour of the integration to bypass the singularity on the left. With the 
e-regularization this means that the contour must pass to the left of —e. As we have seen, 
such a prescription is very natural also when considering the kinetic term. For consistency 
we should therefore adopt the same prescription for computing the coefficients. This will 
lead to a well-defined phantom term, which we shall study later in Sec. 14.21 Note however, 
that this prescription is possible only for integer r, when there is no branch cut extended 
to the origin. The anomaly in the equations of motion for the 'fractional' D-branes is thus 
inevitable. Alternative strategy for canceling the anomaly would be to consider functions 
G with an infinite number of zeros or poles along the imaginary axis. We will postpone 
such an approach to the future work. 

4 Level expansion 

4.1 General arguments for convergence 

As we have already explained, higher multiple D-brane solutions are associated to func- 
tions G{K) that have a pole singularity at the origin, and therefore it is of utmost impor- 
tance to carefully examine whether such solutions make sense. Up to date, inverse powers 
of K have appeared in the literature in attempts to write formally the tachyon vacuum 
as a pure gauge configuration. It has been argued that such string fields are not to be 
considered as good string fields for the purpose of a gauge transformation and hence the 
non-triviality of the tachyon solution remained unaffected. 

The reason why objects like 1/K are dangerous is very simple: at least formally we 
have 

dae-"^, (4.1) 
and regarded as a string field this is linearly divergent since 

lim e-°^ = |oo) = e~5^-2+^^-4+-|o) (4.2) 

is the well-known sliver state [19]. We thus have to be very careful whenever inverse 
powers of K appear. Let us study in more detail string field of the form 

cK^'BcK"', (4.3) 

where n and m are integers, possibly negative. The coefficients of such a field in the stan- 
dard Fock space basis can be easily computed by calculating the overlaps with Fock states, 
or more systematically using the tools developed in [201 E] ^^^^ reviewed in Appendix [XI 
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The starting point is to compute the coefficients in 

e-"^cfie-^^ce-^^ (4.4) 

and obtain the coefficients of cK^BcK"^ by differentiation with respect to the parameters 
/3 and 7 and setting them to zero, or by integration, depending on the sign of n and m. 

The coefficients are given for reader's convenience in a compact form in Appendix |Al 
To write them exphcitly, one has to choose a basis of states. Most convenient one for 
work with analytic universal solutions in Lg-level expansion is the basis formed by b and 
c ghosts, and total Virasoro generators L acting on the SL{2,Wj invariant vacuum |0). 

It turns that potentially most divergent terms are those that include only ghosts, so 
let us focus our discussion only on states of the form c_n|0), ignoring all other for the 
moment. Instead of listing the coefficients tn of c_„|0) explicitly, let us write down the 
asymptotic behavior for large 7 with fixed /3 (setting a = 0)0 





IT 

8 ~ 


7r3 

9672 ^ 


to = 


1 

2 ~ 


7r2 

1272 + 


tl = 


2 

TT 


vr ( 
272 ^ 


t2 = 


8 


8 



4873 

(1 - 3/32 - 2/33)7r2 



+ 0(7- 



37 



6^3 

+ 0(7 



(3 - 6/32 - 4/33)vr , _4, 



3(1 + 2/3)73 



+ 0(7-^). (4.5) 



It is important to notice that all coefficients up to order 7~2 are independent of /3, and 
therefore they would vanish when a derivative with respect to /3 is taken. Starting from 
^2 the coefficient of 7"^ will develop a nontrivial denominator thanks to the presence of 
uncanceled tangents in the expansion of the c-ghost. 

Inspecting closer the coefficients of Ci|0), Co|0), and c_i|0), we find that in general they 
are unambiguously defined in the string field cK^BcK^ whenever m + n > 1. The reason 
being that as long as one of the integers is positive, the respective K turns effectively 
into a derivative and reduces thus the divergence in the integration over the wedge-width 
variable associated with the other power, if it is negative. Note that the criterion m+n > 1 
can be nicely restated in terms of the L~ weight, which is related to the Co weight. The 
nice thing about the L" weight, is that it is strictly additive in the K, B, c algebra. The 
fields K and B each carry weight one, while the field c has weight minus one. So we can 



^^The expansion for large /3 and constant 7 is needed for the discussion of 'ghost branes'. The expansion 
starts at order l3~^ with coefficients which are polynomial or rational functions of 7. 
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say that the coefficients of ci|0), co|0), and c_i|0) in cK"'BcK"^ are finite as long as the 
total L~ weight is greater or equal to zero0 

On the other hand, coefficients of higher level states such as c_2|0) or c_3|0) might 
be divergent even when the condition m + n > 1 is satisfied. The reason is that for 
these coefficients repeated derivative with respect to /3 does not improve the large 7 
behavior. Therefore these coefficients are divergent for m < —3 regardless of the value of 
n. Similarly, they are divergent for n < —3 regardless of the value of m. 

The divergences can be alternatively seen by computing overlaps with the Fock states. 
For the potentially most divergent piece cK"^^^BcK~^ in the multibrane solution we find 

= / dt- -Ti(e-^tcdcVe-^^cK''~^^Bce-^] (4.7) 

Jo (n-1)! V / 



where in the second line we used as an example = cdcV with V being a purely matter 
field of dimension hy > 0. The c ghost OPE's give us a factor of 1/t^ so that the 
integral is nicely convergent at large t. (There is no issue at small t.) Choosing however 
=: cdcd^cb : we get the integral of the form dtt"'~'^. Therefore the coefficient of 
c_2|0) is divergent for n > 3, in agreement with our previous analysis. Similarly we would 
find divergence for n < — 4. 



4.2 Towards the phantom 

As we have just seen, the coefficients of cK^^^BcK~^ in the Fock state expansion are 
finite for —3 < n < 2, and such a statement remains true even if we replace the powers 
of K with arbitrary functions of K, as long as we do not change the behavior near 
K = So we could just compute the coefficients of Fi{K)cF2{K) BcF3{K) by using the 
Laplace representation of the functions Fi{K) and the knowledge of the coefficients in 
e''^^ cBe~^^ ce~"'^ . Throughout the paper we were faced with analogous computations, 
for which we used the so called s-z trick number of times. We have seen on the examples 
of the kinetic term energy (for more than two branes), EUwood invariant, or the anomaly 
in the equation of motion, that we get the correct answer, if and only if our 2;-integration 
contour running along the imaginary axis bypasses the singularity at the origin from the 
left. On the other hand, straightforward integration over a, /3 and 7 gives the same answer, 
as if the z contour bypassed the origin on the right! For numerical computations, presented 

-'^^There are four exceptions to this rule, the states cKBcK^^ , cKBcK^^, cK^^Bc and cK^^Bc with 
L~ weights —1, —2, —3 and —2 respectively, have well defined absolutely convergent coefficients of ci|0), 
co|0), and c-i|0). 
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in the next subsection, it is more convenient to find the coefficients by integrating over 
P and 7 (for this class of solutions a is simply set to zero), than by the computing the 
integrals over s and z. Therefore, we will now describe how to go around the singularity 
on the left, in terms of the /3,7 integrals. 

General coefficient of cF2{K)BcF^{K) can be computed as 



dp / t/7/2(/3)/(/3,7)/3(7), (4.8) 
Jo 

where /2(/3) and 73(7) are inverse Laplace transforms of F2{z) and F^{z). The function 
fiPil) gives a particular chosen coefficient in cBe~^^ce~^^ . This expression can be 
rewritten in full generality as 

1^']. i^/\F2{z)f{-d.,s+dS)F,{z). (4.9) 

For the coefficients of Ci|0), Co|0), and c_i|0) life is simpler once again, the function 
fi—dz,s +dz) depends on the derivative only through a linear factor, or through a 

simple exponential, such as e^i^^"', which acts as a translation operator. For higher 
level coefficients such as c_2|0) or c_3|0) we get infinite number of exponentials (upon 
expansion of tangents) and this presents extra challenges discussed later. 

In the formula (14. 9 p the line integral along the imaginary axis passes through a sin- 
gularity at the origin (due to the presence of either F2 or F3), and one has to specify a 
prescription. The prescription which gives the same result as (14. 8 p can be found in the 
following way. As we have shown in the previous subsection, the integral (14.80 is abso- 
lutely convergent for — 3 < n < 2, and for simplicity we shall restrict our attention to 
this range only. We can thus multiply the integrand by e~^^'^'^^^. The limit £ — )■ is well 
defined and gives back our integral. Now rewriting this e-dependent integral using the 
s-z trick, we get back (14.90 . but with shifted arguments of the F functions, i.e. Fi{z + e) 
instead of Fi{z). Since the contour runs along the imaginary axis and the singularity has 
been shifted to —e, we see that the contour bypasses it on the right. This remains so, 
even in the limit of vanishing e. 

Our proposal for the phantom term is to exclude or subtract the contribution from 
the origin. Let us demonstrate on a simple example of a regular string field 

^ -Bc-J— (4.10) 



K + a K + b 

how we can easily separate the contributions from residues in —a and —b. Let us assume 
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for definiteness a > 6 > 0. For the contribution from —b we need 

F2{z)f{-d.,s+d.) = / e-^72(/5)/(/3,s-/5)| 



(i/3e-^('^-*)/(/3,s-/3). (4.11) 



For the contribution from —a we need 

f{s+d,,-d.)Fs{z) = - I d7e^V3(-7)/(5 + 7,-7)| 



° rf7e-^(''-'^)/(s-7,7). (4.12) 

-oo 

In the first hne we have used the Laplace representation of function vahd in the domain 
Re2; < —b, and in the second hne we just renamed the integration variable from 7 to —7. 
The whole coefficient with contributions from both residues is then simply 

I poo poo pO 

dse-'" d(5e-^^''~^^f{(3,s- (3)- dse'"' d-f e-^^^-"^ /{s - -f,-f). (4.13) 

Jo Jo J-00 

Upon a change of variable 7 = s — /3, this can be written as 

/ poo 

dse-^'i -/ c//3e-^('^-^)/(/3,s-/3), (4.14) 



where the first integral in the bracket gives the contribution from —b while the second 
from —a. The sum of the two is 

dse-'' dPe-^^^-'^f{P,s-P)= ds dPh{P)Ms-P)f{P,s-P), (4.15) 



which is nothing but (14. 8p . 

So we propose a simple prescription for the phantom. From the naive solution whose 
coefficients are computed using the (3 and 7 integrals, one should subtract a finite phantom 
term which is defined as the contribution from the pole of F3 in the case of the double 
brane. Explicitly it takes the form 

poo poo 

Phantom = / ds d(3 f2{(3)f{(3^s - I3)h{s - (3). (4.16) 
Jo Jo 

To conclude, a general coefficient for the double brane solution is given by 

o poo poo poo 

df3 c?7/2(/3)/(/3,7)/3(7)- / ds d^ h{P)f{P,s - P)h{s - (4.17) 
^0 Jo Jo 
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One of the potential subtleties of this phantom, is that for certain higher level coefficients, 
those that contain higher powers of 

tan ( ^ ) = tan ( ^^-^) (4.18) 

the integrand becomes singular for infinite number of values of (3 = s+| + (fc + l)s. Note 
that none of these values lies in the interval (0, s), which is the range of the naive term 
(14. 8 p and which is therefore fine. When we encounter such singularities, we adopt principal 
value prescription. Numerically we implement it by taking the /3 integral originally over 
the positive real axis into the complex /3 plane, slightly above or below the axis and taking 
the arithmetic mean of these two results. This also cancels the unwanted imaginary parts. 

4.3 Numerical results 

So far, most of our discussion of multiple D-brane solutions in previous sections was rather 
general, in terms of a function G{K), which led us to impose some important conditions. 
Let us now look in more detail at the solution 

^ = c- —Be 1 - ^——^ (4.19) 

for which 

G{z)=(^^J. (4.20) 

As special cases, these solutions reduce to the perturbative vacuum for = and to the 
tachyon vacuum [S] for n = —1. If the tachyon vacuum were a large gauge transformation 
of the perturbative vacuum (which it is probably not in any good sense), this whole family 
could be viewed as given by powers of the tachyon gauge transformation. 

It is quite straightforward to compute any coefficient in the expansion of \E'. For 
example for n = 1, which is a candidate for the double brane, we find the tachyon 
coefficient of Ci|0) to be 

t = 0.372994 - 0.588638 = -0.215644. 

The first term is the naive contribution, the second term subtracts the residue at zero. 

In the following tables we list coefficients from Ci|0) to c_5|0) for values of n between 
n = 2 (triple brane) to n = —3 (ghost double brane). The first table includes only the 
naive or regular term. 
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n 


Cl 


Co 


C-l 


C-2 


C-3 


C-4 


C-5 


2 


0.5373 


0.77577 


0.48429 


0.03048 


-0.70900 


-1.7379 


-3.20121 


1 


0.37299 


0.44803 


0.45303 


0.47698 


0.49104 


0.51026 


0.51867 


-1 


0.28439 


0.24903 


0.24452 


0.25257 


0.26894 


0.29239 


0.32283 


-2 


0.63671 


0.23069 


0.16044 


-0.03167 


-0.13167 


-0.33494 


-0.50203 


-3 


1.37982 


0.09646 


0.41052 


-0.31435 


-0.13341 


-0.75982 


-0.68690 



The second gives the final coefficient with the phantom subtracted. 



n 


Cl 


Co 


C-l 


C-2 


C-3 


C-4 


C-5 


2 


-0.43481 


-0.18127 


0.02353 


-0.32736 


-1.20906 


1.45109 


4.13767 


1 


-0.21564 


-0.06671 


0.00911 


-0.08671 


-0.49105 


0.25643 


1.88705 


-1 


0.28439 


0.24903 


0.24452 


0.25257 


0.26894 


0.29239 


0.32283 


-2 


0.64199 


0.48567 


0.40056 


0.35881 


0.70399 


0.93730 


-0.91004 


-3 


0.97760 


0.64625 


0.48330 


0.30180 


0.92385 


2.27210 


-1.31707 



How can we tell whether the results are meaningful? First thing which may worry the 
reader is that the coefficients do not decay rapidly, in fact, they do not decay at all at 
higher levels. But this is a well-known feature of the non-twist invariant solutions, such 
as the ES solution [D] given on the third line. Our solution f l4.19p contains c on the left, 
and therefore it must be annihilated by the operator c(l). 

Another very positive aspect we can read from the table relates to the EUwood invari- 
ant. The Ellwood invariant measures the change in the closed string one point function 
and it can be easily computed using the conservation laws found in [211 ES]. The laws 
allow one to reduce the computation of ( I\ccV{i)\(l) ) for arbitrary in the universal sector 
to the case of |0) = Ci|0). For a string field of the form 

tlCilO) + t2fe-2CoCi|0) + tsC-llO) + t4^-2Cl|0) + t56-2C-iCo 1 0) + t6&-2C-2Ci |0) + t7&-4CoCi|0) 
+t8C-3|0) + tgL_2b-2C0Cl\0) + tioL_2C-i|0) + tllL_4Ci|0) + ti2L_2^_2Cl 1 0) • ■ ■ (4.21) 

the Ellwood invariant, normalized so that for the tachyon vacuum it gives +1, is given by 

I ih +h + te + 4tg - 3tio - 4ti2 + ■■■)• (4.22) 

The odd level fields (with even eigenvalue of Lq) do not contribute, but many other terms 
do not contribute as well. This formula explains why for the tachyon vacuum solutions, 
known to be behaving well in level truncation, the tachyon coefficients are close to 2/7r. 
Incidentally, in the Cq level truncation, the tachyon coefficient of the Bq solution [3] has 
exactly this value. 

Looking at the second table, we see that the Ci coefficient depends on n roughly in a 
linear manner, and that for each conjectured n + 1 brane solution, the contribution from 
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this coefficient gives between 34% to 51% of the expected Ellwood invariant. This is quite 
encouraging, especially in comparison to the tachyon vacuum solution, for which the Ci 
coefficient accounts for about 45% of the expected value, and we know from other studies 
that it behaves well in many respects. 

For more detailed analysis we have selected the double brane solution. From the Fock 
state coefficients up to level 16 we have found the kinetic term and the Ellwood invariant. 
We give our results in the form of a polynomial in 2;, in which the variable z acts as a level 
counting variable. For the kinetic term, normalized so that the expected value would be 
+1 at 2; = 1 we found 

E = -0.152973 - 0.02926592 + 0.260087^^ + 0.0787223^^ - 1.84353z^ + 0.264623z^ 
-1.25582^ + 2.46205;z^ - 26.0792^^ + 14.8247;z^ - 83.7969z^° + 81.068z" 
-450.531z^2 + 387.393^^^ - 1731. 82^^ + 1673.6^^^ - 6957.27z^^ + ■ ■ ■ . (4.23) 

Following [31 |9] we one may attempt to resum this apparently divergent series using 
Fade approximants. We will not present these results in detail, since we did not find 
them illuminating enough, and we think there is more to understand, than what we have 
been able so far. Anyway, our observations from the Fade analysis are that things looks 
promising at lower levels. The energy has to start negative (i.e. with the wrong sign), 
on general grounds, so is reassuring to see that the value —0.15 is not too far from zero. 
Also the second correction which has to be negative is rather small. Then come finally 
two values which give positive contribution. Up to this level, everything looks rather nice. 
But then the positive contributions which come mostly at even levels (odd powers of z) 
seem to be smaller than the negative contributions from both neighboring levels. This 
seems to imply the wrong sign of the resummed answer, but the change in the expected 
value from level to level is too high to draw any definite conclusion. 

Similarly,for the Ellwood invariant, normalized so that the expected value for the 
double brane is — 2/7r we have obtained 

E = -0.215634 + 0.0091296^2 + 0.184435^^* + 0.205632^^ - 1.23245^^ + 1.38278^^° 
-1.7683z^2 + 5.19292^^^ - 11.6944^^*^ + ■ ■ ■ . (4.24) 

Again, the Fade analysis does not provide much help, so we do not present the details 
here. Optically it seems to provide about one third of the expected answer (with the 
correct sign), but the variability is again very high. For comparison, let us show the 
analogous polynomial computed to the same level, in the same normalization, for the 
asymmetric simple solution of [9] 

Etv = 0.284394 + 0.244516^^ - 0.0375272z^ + 0.0504927/ + 0.0925147/ - 0.0026853^^° 
-0.105459^^2 + 0.00828519^^^ + 0.148443^^^ ■ ■ • . (4.25) 
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In this case the convergence of the Pade approximants at z = 1 to the expected value 
is quite apparent, and up to this level one finds results within about 7% of the expected 
answer. 

What should we conclude from these numerical results? One possibility is that we 
have not gone to high enough levels to see the convergence, or that we have not computed 
our numerical coefficients accurately enough. To this level we had to compute about 1316 
coefficients, each given by quite a slowly convergent double integral. Another possibility 
is that our prescription for computing certain coefficients, such as the one of c_3|0) for 
which the phantom naively diverges, is incorrect. Third possibility is that because of 
specific analytic structure of the energy, as a function of the level counting parameter 2, 
the series is simply not Pade summable, and that in the tachyon vacuum case where it 
did work, we were just lucky. More serious possibility, of course, is that we have identified 
the phantom incorrectly. 

5 Discussion 

In this work we have studied Okawa's family of string field theory solutions depending 
on a single analytic function. We have shown how to compute their energy and the 
closed string one-point functions, and that these two computations agree. Our results 
immediately suggest that functions with a pole at the origin should be interpreted as 
describing multiple D-branes. 

One question one could ask is whether our solutions favor multiple D-branes over 
unphysical configurations with negative number of D-branes. They do, but only a little 
bit. For negative n, the solution (14.191) contains a piece with negative weight. Such states 
are generally problematic in correlators and in the overlaps with Fock states. But we 
have not seen anything wrong with the one and two ghost brane solutions. We expect 
that a more detailed analysis would rule out such solutions. Another question is related 
to the existence of fractional number of D-branes. We certainly do not expect their 
existence in the universal sector of bosonic string field theory. Such solutions are plagued 
by irreparable anomalies as we have shown in Sec. [3l 

We have shown that our proposed multibrane solutions have to come with a special 
prescription, of bypassing a singularity at the origin from the left. We have interpreted 
such a prescription as a sort of finite phantom term, but this is not a unique possibility. In 
Sec. |3]we have seen another type of phantom which arises when we systematically replace 
all K^s with K + e. The latter phantom vanishes in level truncation, but at the same time 
it provides the sole contribution to the EUwood invariant. We have tested the proposed 
double brane solution numerically, but the results we got are inconclusive. Settling this 
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issue is left for a future work. It would be very nice if more regular solutions, which would 
behave well in level truncation, were found. 
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A Coefficients in level expansion 

In this short appendix we would like to remind the reader some results from [201 E] and 
show how to efficiently compute coefficients in the level expansion of 

e-"^c5e-''^ce-^^. (A.l) 

This can be written ClS db state in the Hilbert space 

-Ua+B+.+ l [gf- -" + ^ + ^ 1 +5f- -"-^ + ^ ) - iBcf^ -« + /3 + 7 \^(^ -a-/3 + 7 \1 jg) 
TT "+^+7+^ L \2a + t3 + -y + lJ \2 a + 13 + -y+lj w V2a + /3 + 7 + l/ \2 a + l3 + -y + 1 J j ' " 

where Ur = U*Ur, B = Bq + B^, the star denotes BPZ conjugation, and the remaining 
symbols follow the notation of [3]. In particular 

c{x) = cos(x)^c(tana;), (A. 3) 

Ur = 0^ = 0^ e^^^^e"*^* . . . , (A.4) 



27 



where Un are constants given in [3] . More convenient form of the string field flA.ip for the 
purposes of level expansion is a 'normal ordered' form 



vr 



, l.^f TT — a + /? + 7 \ , l,_/7r — a — /3 + 7 



|0) 



a + /3 + 7 + 1 



^/vr — a + /3 + 7 \~/7r — a — /3 + 7 



^:+;3+7+i^o5(^2a + /3 + 7 + l 



2 a + /3 + 7 + 



- |0). (A.5) 



For example the coefficient of ci|0) can be easily read off 



a + iS + 'y + l 



2t 

(a + /3 + 7 + l)^ 
27r2 



7 + - COS ^l^~r I 



+ a+- COS ^ ^ 

^ 2^ V2a + /3 + 7 + l 



(A.6) 



tan(^.-" + ^ + '^ 



2a + /3 + 7 + l 



tan -- 



2 a + /3 + 7 + l 



' + 7 A \ COS^ -" + /^ + ^ ^ COS^ /3 + 7 



2 a + /3 + 7 + l 



2 a + /3 + 7 + l 



B What is behind the s-z trick? 



Undoubtedly much of our discussion in this paper relied on the not so transparent s-z 
trick. In this appendix we will try to clarify it a bit in a simpler setting. Instead of looking 
at complicated strings fields with ghost insertions, let us apply the s-z trick to a simple 
product of functions. Following exactly the same steps as before we get 



i=l 

In particular for single function we get 



ds 



1=1 



w 



(B.7) 



ds 



dw 



e''"F(z + w) 



(B.8) 



What are the conditions of validity of such an expression? Let us look at an instructive 
example of 

The integral along the imaginary axis can be made into a closed contour by adding a 
non-contributing arch at infinity in Rew < half-plane. This contour integral can then 
be nonzero only if it encircles the pole at w = — (z + 1). Therefore the right hand side 
equals to 



F^z) 



1 



1 + z 



^(Re(z+ 1)) 



(B.IO) 



where 6 is the usual Heaviside step function. This agrees with the original function F{z) 
clearly only in the half plane Rez > —1. For more general rational functions F{z) the 
corresponding Fs.z{z) can be defined analogously using the partial fraction decomposition. 
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The function Fs-z{z) is clearly not holomorphic (although it can be analytically con- 
tinued), but this appears to be more of a virtue in cases where the argument is K and we 
have to compute correlators. 
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